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Dynamics  of  a  Gravity  Gradient  Anchored  Tethered  Space  Antenna 

Introduction 

The  intrinsic  attitude  stability  of  a  vehicle  in  orbit  is 
fundamental  to  the  viability  and  efficacy  of  any  configurational 
concept.  While  virtually  any  spacecraft  can  be  stabilized  with 
sufficient  effort,  designers'  attention  is  most  naturally  drawn  to 
those  configurations  which  require  the  least  such  efforts.  This  report 
begins  an  analysis  of  the  unassisted  attitude  stability  of  a  somewhat 
unconventional  spacecraft  configuration,  but  one  for  which  a  number  of 
practical  applications  have  been  suggested.  The  concept  studied  is  that 
of  an  orbiting  antenna  structure  of  laminal  geometry,  which  flies  with 
the  large  dimension  collector  surface  oriented  to  the  NADIR  direction 
and  is  anchored  in  that  attitude  via  tethered  endmasses. 

Figure-1  illustrates  the  notional  system.  In  addition  to  the 
antenna  body,  two  tethered  endmasses  are  considered  which  can  be 
attached  to  the  antenna  by  any  number  of  tethers.  The  tethers 
themselves  are  considered  massless  tensile  members  with  viscoelastic 
constitutive  character.  The  nonlinear  equations  of  motion  are  derived 
for  the  system  subject  only  to  the  gravitational  attraction  of  an 
inert ially  spherical  primary.  An  equilibrium  solution  corresponding  to 
®  circular  orbit  is  determined  and  the  motion  equations  are  linearized 
about  that  NADIR  oriented  steady  state.  The  linearized  equations  are 
shown  to  be  stationary  and  are  brought  to  a  canonical  second  order  form 
to  which  the  theorems  of  linear  stability  theory  can  be  applied. 


Figure-1.  Gravity  Gradient  Anchored  Tethered  Space  Antenna. 
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System  Definition 


The  idealized  system  S  is  comprised  of  a  rigid  body  & 
and  two  point  masses  Px  ,  (i=l,2).  Particle  Px  is  assumed  to  be  connected 
to  P  by  N-£  tethers  attached  at  arbitrary  positions.  P±  can  be 
connected  to  P2  by  only  a  single  tether.  The  tethers  are  idealized  as 
massless  structural  members  capable  of  exerting  only  tensile  force  on 
the  respective  endbodies  directed  along  the  straight  line  connecting 
their  attachment  points.  The  tethers  are  considered  to  be  of  fixed 
unstrained  length  (i.e.  deployment/retraction  is  not  allowed)  but  they 
can  stretch  longitudinally.  The  constitutive  character  of  the  tethers 
is  taken  to  be  linear  viscoelastic;  essentially  they  behave  as  massless 
parallel  spring-damper  pairs.  The  system  is  taken  to  be  free  in  space 
subject  only  to  the  gravitational  attraction  of  a  large  inertially 
spherical  primary.  The  assembly  possesses  twelve  degrees  of  freedom  and 
constitutes  a  flexible  multibody  system,  albeit  a  relatively  simple 
one.  Figure-2  shows  the  idealized  system. 


Figure-2.  Flexible  Multibody  System. 

An  inertial  reference  frame  denoted  9i  is  established  with  its 
orijin,  point-I,  at  the  centroid  of  the  gravitational  primary.  The 
position  vector  Rc  locates  the  centroid  of  S ,  point-c,  relative  to 
point-I.  An  orbital  reference  frame,  ?0,  is  defined  having  its  origin 
at  point-c  and  orthogonal  unit  basis  vectors  ox,  o2.  o3  .  An  additional 
reference  frame  9r  is  defined  as  fixed  to  *  ,  with  its  origin  at  the 
mass  center  of  &  ,  point-0,  and  unit  basis  vectors  i,  j,  k .  Figure-3 
illustrates  S  and  the  basic  reference  points,  reference  frames  and 
vector  geometry.  The  position  r  locates  point-0  relative  to  point-c. 
The  position  of  P±  relative  to  point- 0  is  given  by  the  vector  p.j,  . 
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Figure-3.  Reference  Points,  Frames  and  Position  Vectors. 

Denoting  the  masses  of  ^  and  as  mR  and  m-^  respectively,  the 
total  mass  of  S  is 

2 

mT  =  mR  +  J]  mi 
i  =  l 

Definition  of  the  center  of  mass  of  S  provides  the  relation 

2 

mTJ  +  X  miPi  =  0  (1) 

i  =  1 


System  Momenta 

Denoting  the  position  vector  from  point- 0  to  an  arbitrary- 
material  point  of  ?  as  p,  the  linear  momentum  of  %  is 

-r  rd-_„  d 

Lr  =  I  — — (  Rc+r  +  p)  dm  =  mR  —  (  Rc  +  r ) 
i  dt  dt 


(2) 


The  linear  momentum  for  JR^  is  simply 


Li  =  mi  —  (  Rc  +  r  +  Pi  ) 
dt 


The  total  linear  momentum  of  the  system  is 

2 

L  =  LR  4-  ^  Lj  =  mT  Rc 
i=l 


(3) 


(4) 


where  the  notation  v  denotes  the  time  derivative  of  the  vector  v  as 
observed  in  the  inertial  frame. 
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The  moment  of  momentum  of  S  with  respect  to  point -c  is  defined 
as 


H=  J(r  +  p)x(Rc  +  r  +  p)  dm  +  ^  mi(r +  pi)x(Rc  +  f +  pi  )  (5) 

*  i=l 

Introducing  the  angular  velocity  of  relative  to  »  ©r  ,  we  have 
from  the  law  of  Coriolis 


P  =  ®R*P  (6) 

Expanding  Eq. (5)  and  using  Eqs. (1)  and  (6)  we  obtain 

2  2 

H  =  Ir  •  ®r  +  2  Hij  Pi  *Pj  (7) 

i=l  j=l 

where  we  have  introduced  the  centroidal  inertia  dyadic  of  X  , 

Ir  =  J(p2I-pp)dm 

X 

with  1  being  the  unit  dyadic,  and 


c  mi 

Hij  -  m±(  Sij— —  )  = 

IiIt 

6±j  is  the  Kronecker  delta.  The  law  of  Coriolis  allows  us  to  write 

^  3*  d  _ 

Pi  -  —  Pi  +  Or  x  Pi  =  Pi  +  ©R  x  p±  (8) 

where  the  notation  of  a  circle  over  a  vector  implies  the  time 
derivative  of  that  vector  observed  in  the  frame  7r  .  Substituting  Eq. (8) 
into  (7),  expanding  and  collecting  terms,  the  systems  centroidal  moment 
of  momentum  can  be  written  as 

„  _  2  2 

H  =  le  •  SR  +  j  Pi  x  Pj  (9) 

i  =  l  j=l 

where 

2  2 

i©  =  Ir  +  Y,  2  ^i  j  (  Pi  'Pj  i  -  Pi  Pj  ) 
i  =  l  j  =  l 
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is  identified  as  the  centroidal  inertia  dyadic  of  the  system. 


System  Dynamics 

The  forces  acting  on  the  bodies  of  the  system  arise  from  the 
tether  interconnections  and  the  gravitational  attraction  of  the 
primary.  It  will  be  necessary  to  develop  explicit  expressions  for  these 
forces,  but  presently  we  shall  only  define  them  symbolically.  Let  f±  be 
the  gravitational  force  acting  on  P± ,  and  let  fR  be  the  resultant  of 
the  gravitational  force  distribution  acting  on  ^ .  We  define  g  to  be 
the  resultant  moment  about  point-c  of  all  gravitation  forces  acting  on 
S.  Acting  internal  to  the  system  are  the  tether  forces.  Define  4  to 
be  the  resultant  of  all  tether  forces  acting  on  P±  .  Since  the  particles 
connect  to  each  other  as  well  as  to  *  we  distinguish  the  components  of 
4  as  ; 


4  -  Tpi  +  TRi  (10) 

where  TPi  is  the  force  on  P±  due  to  its  connection  to  Pj ,  (i#j)  and 

TRi  is  the  resultant  force  on  p±  from  all  the  tethers  connecting  it  to 
X  .  With  these  definitions  we  can  write  the  dynamic  equilibrium 
equations . 

Newton's  second  law  applied  to  P requires 

mi(  Rc  +  ^  +  Pi  )  =  4  +4 

Using  Eq. (l)  to  eliminate  r ,  the  above  equation  becomes 

2 

miRc  ^  j  Pj  =  4  +  4 
3=1 

Expressing  the  time  derivatives  of  pj  in  terms  of  those  observed  in  ^ 
we  write 


miRc 


-  c 


^  OO 

:±x©R  +  £  n±j  pj 

j=l 


=  4+4-14 


(ii) 


where 


o 

Ni  =  2®RXci  -©RXt^xci) 


(12) 
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and 


2 

5i  = 

3=1 

Newton's  second  law  applied  to  the  system  as  a  whole  states  that 

„  2 

mTRC  =  %  +  5^4  (13) 

i=l 

The  balance  law  for  the  moment  of  momentum  of  S  expressed  in 
terms  of  the  system  centroid  states 


H  =  H  +  ©R  x  H  =  g 

Substituting  Eg. (9)  into  the  above  equation  and  carrying  out  the 
indicated  differentiation,  we  obtain 

-A  „ 

I©  '  ®R  +  2-t  Ci  x  pi  =  §  ~  Ns  (14) 

i=l 

where 

O 

Ns  =  I©  •  ©R  +  ©n  x  H  (15) 

Equations  (11) ,  (13)  and  (14)  constitute  four  independent  vector-dyadic 
motion  equations  for  the  system. 


Gravitational  Forces  and  Moments 

In  this  section  we  develop  approximate  expressions  for  the 
gravitational  forces  acting  on  the  system  components  as  well  as  the 
resultant  moment  of  those  forces  about  the  system  centroid. 

The  gravitational  force  acting  on  is  given  by 

f.  =  gig*  +  £  +  Pi  ) 

I  Rc  +  r  +  pi  | 3 


where  p  is  the  gravitational  constant  of  the  primary.  The  resultant  of 
the  gravitational  force  distribution  acting  on  &  is 


f gfc+r+pl 
*  |RC  +  r  +  p|3 


(17) 
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Consider  the  expression 


|Rc+A|_3=  [  (Rc+A)-(Rc+A)]”^  =  R“3  (  1  +  2  )-% 

Rc  Rl 

For  the  problems  under  consideration  here  Rc  » A  .  Applying  the  binomial 
expansion  and  retaining  terms  to  only  first  order  in  (A/Rc),  we 
write 


|RC+A|  3s  rc3(1-35s^.)  (18) 

Rc 


Using  the  approximation  of  Eq. (18)  in  Eqs. (16)  and  (17)  we  obtain 


4  s  {  °3  “  t  (r+Pi)  -  363  63-(r+pi)] }  (19) 

Rc  Kc 


?r  =  —4^-  [  o3  -  -i-  (  r  -  3o3 63  -r  )  ]  (20) 

Rc  Rc 


where  the  unit  vector  63  is  established  in  the  NADIR  direction; 

°3  =  -Rc  /  Rc  *  The  resultant  of  the  moments  of  the  gravitational  forces 
acting  on  S  about  point-c  is 


g  = 


-  J(r+p)x 

* 


p(  Rc  +  r  +  p  ) 
|RC  +r  +  p|3 


2 

dm  +  ^(r+PiJxfi 
i=l 


(21) 


Using  the  approximation  of  Eq. (18)  in  Eq. (21)  and  ignoring  terms  of 
second  or  higher  order  in  (A/Rc)  in  relation  to  terms  of  zero  or  first 
order,  we  obtain 


g  £ 


r  x  03  + 


-y  63  x  JR.o3  +  (r+p±  )  x  f± 

R<=  i=l 


(22) 


where 


5r  -  Ir  +  mR(  r2  1  -  r  r  ) 

is  the  inertia  dyadic  of  X  with  respect  to  point-c.  Substituting 
Eq. (19)  into  (22)  and  simplifying  we  achieve 

_  3\i  * 

<3  =  — r03  XIe  -o3  (23) 
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Tether  Forces 


Figure-4  illustrates  the  tether  forces  acting  within  S .  Consider 
first  the  single  tether  connecting  to  ^  .  Define  sp  to  be  the 
position  vector  from  to 


5P  =  P2  “  Pi 

The  distance  between  the  particles  is 

SP  =  (  sp  ■  Sp  ft 

and  the  unit  vector  in  the  direction  of  Sp  is  denoted  sp  .  The 
unstrained  length  of  the  tether  connecting  the  particles  is  denoted  as 
the  constant  lp .  kp  is  defined  as  the  longitudinal  stiffness  of  the 
tether  and  is  constant  assuming  uniform  material  and  geometric 


<s>, 


Figure-4.  Tether  Forces. 

properties  (kp  =  EA/lp,  where  E  is  the  modulus  of  elasticity  anri  a  is 

the  effective  cross-sectional  area) .  The  inter-particle  tether  forces 
are 


Tpi  =  -Tp2  =  Tp  sp 


(24) 
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where  the  tether  tension  is  defined  as 


kp  ^  SP  ^p  ^  "^p  ^  0 

0  otherwise 


(25) 


The  tension  must  be  non-negative.  The  stretch  rate  coefficient  dp 
governs  the  viscous  damping  of  the  tether. 

Let  tij  designate  the  jth  tether  connecting  p±  to  j?  ,  (i=l,2), 
(]=i|2, . . . ,Nj)  .  tjj  attaches  to  K  at  point- A^j  .  The  position  vector 
aij  locates  Aij  relative  to  0  (see  Figure-4).  The  position  vector  from 
Ai j  to  P^  is  denoted  as  j ,  and  its  associated  magnitude  and  unit 
vector  are  sAj  and  respectively.  Let  T±j  be  the  tension  force 
developed  in  x± j ,  and  define  k±j  and  d±j  as  the  tether's  stiffness  and 
damping  coefficients  respectively.  Then  we  define  the  constitutive  law 
for  Tij  as 


Tij  = 


^ij(sij  lij)  +  Tj^j  ^  0 


otherwise 


(26) 


where  lij  is  the  unstrained  length  of  Tij .  The  force  applied  to  P±  by 
Tij  is 


lij  —  lij  ®ij 


and  the  resultant  of  all  such  forces  is 


iRi  =  ~  Tij  Sij 

j=l 


Hi 


-S’ 

j=l 


( aij-pi ) 


where 


(27) 


sij  “  [  (pi  —  aij )*(pi  —  aij )  ]/2 

Equations  of  Motion 
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Having  developed  explicit  expressions  for  the  internal  and 
external  forces  acting  on  S  we  can  now  complete  the  specification  of 
the  nonlinear  motion  equations  embodied  by  Eqs. (11),  (13),  and  (14). 

Consider  first  the  translational  equilibrium  equation  for  the 
system,  Eq. (13) .  Substituting  Eqs. (19)  and  (20)  into  (13)  we  find  that 


Rc  ~ 


(28) 


Equation  (28)  is  recognized  as  the  classical  two-body  motion  equations. 
To  the  level  of  approximation  made  here,  the  motion  of  the  systems' 
centroid  is  independent  of  the  other  degrees  of  freedom.  Since  the 
solution  of  Eq. (28)  is  fully  known  we  shall  not  consider  it  further. 

The  appearance  of  Rc  in  Eq.  (11)  can  be  eliminated  by  substituting 
Eq. (28)  into  (11) .  Further  simplification  can  be  achieved  by  using 
Eq.  (1)  to  eliminate  the  explicit  appearance  of  r  in  Eq.  (19)  ,  and  in 
turn  from  Eq. (11) ,  leading  to 

oo  ( 

Pj  -  cix©R+  =  Tpi  +  TRi  -  -^-  (  1-36363  )-Ci  -  Ni  (29) 

j=l  Rc 


To  obtain  the  scalar  form  of  the  motion  equations  all  the  vector 
and  dyadic  quantities  appearing  in  Eqs. (14)  and  (29)  are  resolved  into 
their  components  expressed  in  .  The  3x1  column  matrix  of  the  scalar 
components  of  a  vector,  say  pA  ,  resolved  in  9r  are  denoted  by  the 
under-bar  notation,  p±  .  The  inertia  dyadics  I©  and  iR  resolved  into 
components  in  ^  ,  are  represented  by  the  3x3  matrices  [Ie]  and  [IR]  . 

We  also  introduce  the  matrix  cross  product  notation; 


where 


V  X  u 

0 

VXU 

'V1] 

0 

-v3 

v2  " 

v  = 

v2 

yx  = 

v3 

0 

-Vi 

,  v3  , 

.  -v2 

Vi 

0 

With  this  notation  established  the  nonlinear  motion  equations  are 
expressed  as 


[I©  ] 

<2 

_ I 

•  > 
“r 

(  \ 
g 

f  0  > 

rNs> 

mm 

U12U] 

Hi 

= 

2i 

+ 

-PI  +  — R1 

- 

Hi 

P2l(  U 

P22  ( 1  ] 

,22  j 

,  ^P2  +  ^R2  > 

(30a) 
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where  we  have  introduced  the  3x3  identity  matrix  [  1  ] ,  and 

2 

— s  =  ®r  [  xe  ]  ®r  +  2  X  (  -R£i  ~  Si-R E  1  ]  }  Ei 
i=l 

2 

2*i  =  2  ®R  Z  ^ij  Pj  +  ®R  ®R  Si 
j=l 

3p  x 

£  =  — r  °3  C  !©  ]  23 

Rc 

2i  =  -4([U-3o3o|)ci 
2  2 

t  I©  ]  =  [  IR  ]  +  £  Z^ElEjlH-EiSj  ) 
i  =  l  j=l 


(30b) 

(30c) 

(30d) 

(30e) 

(30f) 


Attitude  Kinematics 

A  primary  interest  of  this  study  is  the  attitude  behavior  of  the 
system  relative  to  a  suitable  orbital  reference  frame.  The  orbital 
frame  selected  is  the  usual  one  used  for  spacecraft  attitude  dynamics, 
namely  a  local -vertical  local-horizontal  frame  with  axes  defined  by; 


03  =  - 


RrXR_ 

°2  =  — - ~ 


I  Rc  x  Rc  I 


°1  =  02x°3 


The  motion  of  the  orbital  reference  frame  is  completely  determined 
from  the  solution  of  the  two-body  problem  governing  the  motion  of  the 
system  centroid  Eq. (28) . 

The  angular  velocity  of  ?0  relative  to  Pi  denoted  ©G ,  is  known 
in  terms  of  the  true  anomaly  v(t) 

©o=-vo2  (31) 

The  angular  velocity  ^  relative  to  ?0  is  denoted  as  Q  ,  so  we 

have 


®R  -  ®0  + 

The  angular  acceleration  of  5^  is  noted  as 


(32) 
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(33) 


O 

©R  =  ©0  +  ft  +  ©Q  X  Q 

The  attitude  of  Sk  relative  to  ?Q  is  described  in  terms  of  a  1- 
2-3  Euler  angle  sequence.  The  signs  of  the  single  axis  rotations 
defined  below  are  positive  according  to  the  right  hand  rule.  Let  ?0  be 
an  intermediate  frame  achieved  from  ?0  by  a  rotation  of  angle  0X  about 
the  o1  axis.  Let  frame  be  achieved  from  by  a  rotation  of  angle 
02  about  the  axis  o2  .  5^  is  achieved  from  ^  by  a  rotation  of  angle  03 
about  the  k  axis.  The  direction  cosine  matrix  transforming  vector 
components  from  to  is  obtained  by  concatenating  the  sequence  of 
single  axis  rotations  yielding 


C02  C03 

c0i  s03  +  s03  s02  c03 

:  S03  s03  -  C03  s02  c03 

(cRoJ  ~ 

-C02  s03 

C0!  c03  -  s03  s02  s03 

:  S03  C03  +  C03  S02  S03 

s02 

-s03  C02 

c03  c02 

where  we  have  used  the  shorthand  notations  cO-^cosOi,  s0i=sin0i.  The 
direction  of  the  transformation  of  Eq. (34)  corresponds  to 


/  r  > 
1 

f  - 

°1 

j 

=  (  Cro  l 

°2 

UJ 

<  ®3  j 

The  angular  velocity  of  with  respect  to  ?Q ,  expressed  in 
terms  of  the  Euler  angle  rates  is 


Q  —  0iO2  +  02  ©2  +  63  k 

Resolved  into  components  referred  to  5^  we  have 


002  c03  s03  0 

o  = 

—c02  s03  c03  0 

®2 

S02  o  1 

®3 

Inverting  the  above  equation  yields 


C03/c02 

-  S03/c02 

0 

= 

s03 

C03 

0 

®3 

-  C03t02 

s03 102 

1 

or  more  compactly. 


(35) 
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0  =  [n  ]_1  n 

where  t02  =  tan  02  ,  and  the  singularity  at  02  =  ±  n/2  is  noted. 

Nonlinear  Equations  -  First  Order  Form 

Equations  (32)  and  (33)  can  now  be  expressed  in  scalar  form  as 

®r  =  0  -  v  [  CRo  ]  e2  (36) 

Sr  =  Q  ~  V  [  CRo  ]  e2  -  v  [  CRo  ]  e£  [  CRo  ]T  Q  (37) 


where  e|  =(010). 

Equations  (30a)  are  now  written  as 


[M]{v}  =  {g}  +  {T}  -  {  N  }  (38) 

where  {  v  }T  =  (  p?;  )  -  Here  we  introduce  the  vector  of  generalized 

coordinates 


{q} 


'  0  ^ 


ill 

^2; 


The  mechanical  state  vector  is  defined  as 


m 


(V) 

(q)  . 


The  equations  of  state  for  the  system  can  now  be  expressed  as 


(  Y  ) 


[M]  1(  (g)  +  (T)  -  (N)  ) 

(  V  }  ■ 

tnj_1(  ©R  +  V[CRo]e2 

{ q )  J 

- 

Ei 

E2 

{  f(Y,t)  } 


(39) 


Equations  (28)  and  (39)  complete  the  specification  of  the  nonlinear 
system . 
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Steady  State  Solution  and  Linearization 


We  propose  the  existence  of  a  steady  state  solution  corresponding 
to  the  system  configuration  shown  in  Figure-1  with  the  system  centroid 
in  a  circular  orbit.  The  equilibrium  solution  considered  has 


0(t)  =  0  ,  (constant)  (40a) 

£i(t)  =  p?  ,  (  i  =  1,  2  )  ,  (p°  constant)  (40b) 

Rc  =  -Rc  °3  •  (Rc  constant)  (40c) 

©o  =  “®c  °2  '  (®c  =  v  constant)  (40d) 


For  this  case  the  configuration  is  fixed  relative  to  ?0 ,  and  5^  is 
parallel  to  .  The  explicit  relationship  between  and  need  not 
be  considered  for  this  problem. 

To  study  the  attitude  stability  of  the  system  the  motion 
equations  are  linearized  about  the  steady  state  solution.  Let 


0(t)  =  A0(t) 
Hi(t)  =  P°  +  Ap±(t) 


(41a) 

(41b) 


The  A  quantities  measure  displacements  from  the  equilibrium  state  and 
are  considered  to  be  of  sufficiently  small  magnitude  that  terms  of 
second  or  higher  order  in  them,  or  their  derivatives,  can  be  ignored  in 
relation  to  those  of  first  (or  lower)  order. 

Substituting  Eq. (41a)  into  (34)  yields 


t  CRo  ]  3  [  1  ]  -  A0* 


To  first  order  Eq. (35)  gives 


Q  3  Ae 

Equations  (36)  and  (37)  simplify  to 

— r  =  ~ ®c  §2  +  A0  —  <ac  e^  A0 
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Continuing  our  expansions  of  terms  into  zeroth  and  first  order  terms  we 
record 


with 


Si  =  Si+ 

j=l 


Si  = 

j=l 


(42) 


and 


[I©]  S  [l|]  +  [1^] 


(43) 


where 


2  2 


(  I©  J  -  [  Ir  ]  +  ^  X  l*i j  £  (P°)T  P °  t  1 1  ~  P°  (p°)T  ] 
i=l  j=l 


2  2 


tI©]=  J  £nijn(p°)TAp.  -HP^APi  )  t 1 3  -  <P?<*Pi>T  +  Ap.;<p2>T  )1 

i  =  lj=l 


Expressing  63  in  ^  we  have 


o3  =  e3  +  e*  A0 

where  e|  =  (001).  Using  the  above  results  Eq.  (30d)  is  approximated  as 

9  =  g°  +  g1  (44) 


where 


g°  =  3©*  e3  [  I©  ]  e3 


g1  =  3©|  [  [  l&  ]  -  ( [  i%  ]  e3  )x  ]  e*  A0 


3©^  ^  (  £3  £i  §3  +  §3  c®  §3  ]  Api 


i=l 


Note  that  we  have  used  the  relation  p/R^  =  ©c  • 


For  Eq. (30e)  we  derive. 


(45) 


2i  s  4  +  si 


where 


2°  =  -®|  ( 1 1  ]  -  3  e3  e|  )  c9 

Si  =  3®i[e3(c°)T+e*c°]e*  A0  -  ®|  £  |t±j  ( [  1  )  -  3  §3  e|  )  Ap. 

j=l  3 

To  form  the  linearized  approximations  for  the  inter-particle 
tether  force  we  note  the  following 


§p  s  §p  +  §p 


where 


4 


(£2-e;> 


4  =  <*E2“*El  > 


Also 


s„  S 


4  +  4 


where 


=  m oVTeo  .  _i  _  <4)T4  ,  -i  i  sj, 

V  (-p'  §p  '  sp - 5 —  f  <sp)  =  — - p 


(s 


o  >2 


Expressing  sp  resolved  in  as  u  we  have 


hp  =  4  +  4 


~P  *SP*  -p  ;  4  ~  (sp*  3t  (sp)2[ 1  ]  ~  §p  (gp)T  ]  (  Ap2  -  Apj^  ) 

Using  the  expansions  above  Eq. (24)  is  now  written  as 


Ipx  =  Ipi  +  lli 


where 


(46) 
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Tpi  =  -(  28-q  -  1 )  [  Kp  ]  (  Apj  -  Ap2  )  -  (  28i;L  -  1  )  [  Dp  ]  (  Aj^  -  Ap2  ) 


and 


[Kp]  =  kp[(l-±E-)  (u£)Tu£[l]  +  ^uJ(u£)T]  =  [Kp]T 


[  Dp  ]  =  dp  (u")1  =  [  Dp  ]L 

Development  of  the  linearized  expressions  for  the  tether  forces 
connecting  P^  to  P  requires,  as  above,  the  following  relations 


sl.  ~  JL 

siJ  n 
s 


ID 


Hij  = 


Sij  * 

5ij+Api 

0  0 

;  Sij  =  Pi -Sij 

3  si)  +  s5j 

;  si  j 

=  [§ij(aij-2p°  )  +  (p°)T 

£/*Pi 

;  (s  ±j)_1  = 

1  _  sij  .  Sij 

s°j  (s°j)2  '  s±j 

ii  . 
4i  ' 

-i  = 

r3  [  ( s§  j  )2 1 1  ]  -  2  s9  j  (s9  j  j' 

=  * 
-13  —13 


With  the  above  expressions  it  follows  from  Eq. (27)  that 


— Ri  =  SRi  +  — Ri 


Ni  1-  • 

-Ri  =  ~^Li  ki  j  * 1  *  -i  j  *  $Ri  =  “  C  KRi  ]  “  [  °Ri  ]  AjK 

j=l  sij 


% 

[Kr±  ]  =  2,‘sij>  3  C(sij)2  ki  j  <s?j  —  li  j )  f  1  ]  +  2  k±  j  li;j  s9.  (s9j)T  ]  -  [KRi 


j=l 


t  DRi  ]  =  2  £  (  4  j  r2  di  j  s?  j  (s?  j  )T  =  [  DRi  f 
j=l 


(47) 


Expanding  Eq. (30b)  into  zeroth  and  first  order  terms  we  obtain 


(48) 


where 


Ns  =  +  Ng 


N2  =  ©c  ®2  1 1®  1  e2 


N s  =  ®c 


[  ( [  le  ]  §2  )x  -  §2  [  I®  ]  ]A0  +  2<»c  XtsIc?[l]-e2(cO)T]Ap. 


-  m2  [  ( [  1%  ]  e2  )x  -  e*  [  1%  ]  ]  e*  A0  -  o2c  £  [  e* c<?  ef  +  e*  c?  e*  ]  Ap. 


Expanding  Eq. (30c)  in  a  manner  similar  to  that  above. 


where 


N±  =  N?  +  n| 


NO  rA2  x  x  0 

_j_  —  C0c  §2  §2 


— 1  =  ®c  [  (  ®2£i  )X  +  S2  <S°)X  J  A0  -  2©c  e^  Hij  Ap. 

j  =  l 

“  ®ct(  §2®i  f  +  §2  <£i)X  1  §2  +  ©I  e*  e2  £  Mij  Ap. 

j  =  l 

The  vector  of  generalized  coordinates  is  decomposed  as 

(  q  )  =  (  q°  }  +  {  Aq  } 


where 


(q°)  = 


;  {  Aq  }  =  Ap^ 

Ap, 


The  vector  of  accelerations  appearing  in  Eq. (38)  is  related  to  the 
derivatives  of  the  generalized  coordinates  by 

{  v  }  =  {  Aq  }  +  [ B  ]  {  Aq  ) 


where 
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[B]  = 


-e>c§2 

[0] 

(0] 

[0] 

[0] 

[0] 

[0] 

[0] 

[0] 

Substituting  Eqs.  (42-51)  into  Eq.  (38)  and  retaining  terms  to  first 
order  only,  the  motion  equations  can  be  written  as 

[  M0  ]  {  Aq  }  +  [  M0  ]  [  B  ]  {  Aq  }  =  {  g°  }  +  {  g1  }  +  {  T°  }  +  {  T1  }  -  {  N°  }  -  {  N1  }  (52) 
where 


[Mq  ]  = 


[  I©  ]  (  c°  )x  (  c°  )x 

-(  £1  f  Vll  ( 1  ]  P12  [  1  ] 
-(  >*  H21  [  1  1  P22  [  1  1 


I  g°  1T  =  (  g°T  g if 
T  =  (  oT  if 


(T°  ) 


{  N°  }T  =  £  N f  N  f 


<**-[*'  sf) 

,t1,t = CeT  #  t-) 

(N1)1  =  Nf  Nfj 


All  the  first  order  terms  on  the  RHS  of  Eq. (52)  are  proportional  the 
Aq's  and  their  time  derivatives.  Since  the  steady  state  solution 
prevails  when  all  the  Aq  quantities  are  zero,  the  equilibrium 
conditions  to  be  satisfied  by  the  steady  state  are; 


f 


h2 


\ 


V 


N$ 


4 


J 


o  > 

(  \ 

g 

0 

2? 

rpO 

—1 

O  <M 

w 

mO 

ll2  J 

(53) 


The  first  partition  of  Eq. (53)  requires  that  5^  (parallel  to  ?G) 
correspond  to  the  centroidal  principal  axes  of  the  system.  The  second 
and  third  partitions  of  Eq. (53)  equilibrate  the  steady  state 
accelerations  of  the  endmasses  with  the  net  gravitational  and  internal 
tether  forces. 
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Expanding  the  first  order  terms  on  the  RHS  of  Eq. (52)  and 
recognizing  Eq. (53) ,  the  linearized  motion  equations  can  be  written  as 


t  M0  ]  {  Aq  }+([G]  +  [D]){  Aq  }  +  ([K]  +  [C]+[g]){  Aq  }  =  {  0  } 

where 


t  ([I®]e2)x-[I^]e^  —  tl© )  3  |2[e|c?  [1]-  e2(c£)T] 

2  (§2-2  CD  -  e2(c^)T] 

M  =  ©c 

[  ex  (c°)x  +  (c°)xex  +(ex  c°  )x]  1  -  2pn  ex 

_  2Pi2  ex 

^  <=1  ^  +  i=g  i-sg + £g  :  -2F21?* 

~  2^22  ®2 

[0] 

[0] 

[0] 

[  D  ]  = 

[0] 

C  Dp  ]  +  [  DR1  ]  i 

-  [Dp] 

;°] 

-  [Dp  ] 

[  Dp  ]  +  [  Dr2  ] 

[  g  ]  =  3©2 


-3  I  [  -3Sl§3  +  §3  2i®3  1  ;  t  §3 2° §3  +  §3  ] 


"(.-3.(9?.>T+..-3.-°.).-3.j..3  Pll(  [  1 1  “  ^3*3  )  |  \  Pl2<  [l]-3e3eT> 

®3  (2°  +  §3  9°  )  §3  ;  3  P2l(  [  1 )  “  3f3 ^  )  ;  j  fl22<  [  1  ]  ~  3§3 e|  ) 


[  ]e2  f-  ex[I§,]  ]  ex  !  [  +  e|  c°ex  ] 

[  ^  9°  §2  +  §2  92  ®2  1 

[  C  ]  =  -<o| 

[ex(c0)x+(^c?)x]ex;  -mexex 

“  Ul2®2®2 

[  ®2  <?2  ^  +  <  ®2  -2  **  ^  §2  1  “  P21®*  -2 

-  F22®2®2 

[  0] 

[0] 

[0] 

[0] 

[  Kp  ]  +  [  KR1  ] 

-  [Kp  ] 

_[  °] 

-  [  Kp  ] 

[  Kp  ]  +  [  Kr2  ] 

The  coefficient  matrices  in  Eq. (54)  are  each  constant  and  trace 
their  respective  origins  to  specific  dynamical  effects; 

[  Mq ]  -  generalized  mass  matrix 

[K]  -  elastic  stiffness  matrix  from  the  tethers 
[  D ]  -  viscous  damping  matrix  from  the  tethers 
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t  G ]  -  gyrie  matrix  capturing  Coriolis  effects  from  system  rotation 
[  C  ]  -  centrifugal  matrix  from  system  rotation 
[  g ]  -  gravitational  effects  matrix  for  orbiting  system 

It  is  immediately  evident  that: 

[  Mo  ]  =  [  Mo  f  ;  [  D  ]  =  [  D  ]T  ;  [  K  ]  =  [  K  ]T 
While  less  apparent  it  can  be  shown  that: 

[  G  ]  =  -  [  G  ]T  ;  [  C  ]  =  [  C  ]T  ;  [  g  ]  =  [  g  ]T 

With  these  properties  noted,  Eq. (54)  are  seen  to  be  of  a  second  order 
stationary  damped  gyric  form1.  Important  aspects  of  the  attitude 
stability  of  the  orbiting  assembly  can  be  determined  from  the  further 
study  of  the  properties  of  the  coefficient  matrices. 


1  Hughes,  P.C.,  'Spacecraft  Attitude  Dynamics',  Appendix-A,  John  Wiley  & 
Sons,  1986. 
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